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CIh Multifield models of inflation with nonminimal couplings are in excellent agreement 



< 
(N 



O 

H 



with the recent results from Planck. Across a broad range of couplings and initial 
conditions, such models evolve along an effectively single-field attractor solution and 



/^ predict values of the primordial spectral index and its running, the tensor-to-scalar 

>-' ratio, and non-Gaussianities squarely in the observationally most-favored region. 

^L( Such models also can amplify isocurvature perturbations, which could account for 



the low power recently observed in the CMB power spectrum at low multipoles. 



jrt Future measurements of primordial isocurvature perturbations and the tensor-to- 

scalar ratio may serve to distinguish between the currently viable possibilities. 
> 
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t:;J^ Early-universe inflation remains the leading framework for understanding a variety of 

^f) features of our observable universe [U |2]. Most impressive has been the prediction of pri- 

'' mordial quantum fluctuations that could seed large-scale structure. Recent measurements 

K> of the spectral tilt of primordial (scalar) perturbations, n^, find a decisive departure from 

Cd a scale-invariant spectrum |3l H]. The Planck collaboration's value, Ug = 0.9603 ± 0.0073, 

departs from n^ = 1 by more than 5a. At the same time, observations with Planck constrain 

the ratio of tensor-to-scalar perturbations to r < 0.11 (95% CL), and are consistent with the 

absence of primordial non-Gaussianities, /nl ~ [U [5]. 

The Planck team also observes less power in the angular power spectrum of temperature 
anisotropics in the cosmic microwave background radiation (CMB) at low multipoles, i ~ 
20 — 40, compared to best-fit ACDM cosmology: a 2.5 — 3cr departure on large angular 
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scales, ^ > 5° [6] . Many physical processes might ultimately account for the deviation, but a 
primordial source seems likely given the long length-scales affected. One plausible possibility 
is that the discrepancy arises from the amplification of isocurvature modes during inflation 

In this brief paper we demonstrate that simple, well-motivated multifield models with 
nonminimal couplings [3 |8] match the latest observations particularly well, with no fine- 
tuning. This class of models 

• generically includes potentials that are concave rather than convex at large field values; 

• generically predicts values of r and n^ in the most-favored region of the recent obser- 
vations; 

• generically predicts /nl ~ except for exponentially fine-tuned initial field values; 

• generically predicts ample entropy production at the end of inflation, with an effective 
equation of state Wes ~ [0, 1/3]; 

• generically includes isocurvature perturbations as well as adiabatic perturbations, 
which might account for the low power in the CMB power spectrum at low multi- 
poles. 

We consider this class of models to be well-motivated for several reasons. Realistic mod- 
els of particle physics include multiple scalar fields at high energies. In any such model, 
nonminimal couplings are required for self-consistency, since they arise as renormalization 
counterterms when quantizing scalar fields in curved spacetimes [9]. Moreover, the non- 
minimal coupling constants generically rise with energy under renormalization-group flow 
with no UV fixed-point [10], and hence one expects |^| ^ 1 at inflationary energy scales. 
In such models inflation occurs for field values and energy densities well below the Planck 
scale [71 El E]. Higgs inflation [llj, which must be treated as a multifield model since the 
Goldstone modes remain in the spectrum at high energies in renormalizable gauges [S] , is a 
particularly elegant example of this class of models. 

We demonstrate that models of this class exhibit an attractor behavior: over a wide range 
of couplings and fields' initial conditions, the fields evolve along an effectively single-field 
trajectory for most of inflation. The attractor solution means that for most regions of phase 



space and parameter space, this general class of models yields values of n^, r, the running of 
the spectral index a = dus/dlnk, and /nl in excellent agreement with recent observations. 
Unlike several models with concave potentials analyzed in |1], meanwhile, these models 
should support efficient (p)reheating, readily producing significant entropy at the end of 
infiation. Whereas the attractor behavior creates a large observational degeneracy in the r 
vs. Us plane, the isocurvature spectra from these models depend sensitively upon couplings 
and initial conditions. 

In the Jordan frame, the fields' nonminimal couplings remain explicit in the action. 
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where quantities in the Jordan frame are marked by a tilde. Performing the usual conformal 
transformation, g^„{x) — > g^„{x) = 2M~^ f {(p^ {x)) g^y{x)^ where Mpi = (SttG)^-*^/^ = 2.43 x 
10^^ GeV is the reduced Planck mass, we may write the action in the Einstein frame as [7j 
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The potential in the Einstein frame, 1^(0^), is stretched by the conformal factor compared 
to the Jordan-frame potential: 
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The nonminimal couplings induce a curved field- space manifold in the Einstein frame, with 
a metric given by 
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where fj = df /d(j)^ . For the models of interest, we take 
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Here we consider two-field models, I, J = (p, x- 

As emphasized in [71 El [H], the conformal stretching of the effective potential in the 
Einstein frame, given by Eq. ([3]), generically leads to concave potentials at large field 
values, even for Jordan-frame potentials that are convex. In particular, for a Jordan-frame 
potential of the simple form 
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FIG. 1: Potential in the Einstein frame, V{4>^). Left Xy^ = 0.75 A,^, g = X^, ^^ = 1.2 ^0. Right: 
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^y. In both cases, ^/ ^> 1 and < Xi,g < 1. 



Eqs. (Is]) and (IS]) yield a potential in the Einstein frame that is nearly flat for large field 
values, 
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(no sum on J), as the Jth component of (f)^ becomes arbitrarily large. This basic feature 
leads to "extra" -slow slow-roll evolution of the fields during inflation. If the couplings Xj 
and ^j are not equal to each other, 1^(0^) develops ridges separated by valleys [7]. Inflation 
occurs in the valleys as well as along the ridges, since both are regions of false vacuum with 
V ^ 0. See Fig. 1. 

Constraints on r constrain the energy scale of inflation, H{t^)/M^\ < 3.7 x 10~^ [4j. For 
Higgs inflation, with Xj = g = X^ and ^i = ^(/,, Eqs. ([7| and (10) yield 
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Measurements of the Higgs mass near the electroweak symmetry-breaking scale require A^ ~ 
0.1. Under renormalization-group flow, A,^ will fall to the range < A^ < 0.01 at the 
inflationary energy scale. Eq. ([s]) with A^ = 0.01 requires $,^ > 780 at inflationary energy 
scales, and hence ^<^ ~ (9(10) at low energies [12]. For our general class of models, we 
therefore consider couplings at the inflationary energy scale of order Xj,g ~ C(10~^) and 

e/~ 0(103) m- 



(9) 



Expanding the scalar fields to first order, (f)\x^) = ^p\t) + 6(j)\x'^), we find [71 18] 
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We also expand the spacetiine metric to first order around a spatially fiat Friedmann- 
Robertson- Walker metric. Then the background dynamics are given by [7] 
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(10) 



Vtif' + 3Hip' + G'^^Vk = 0, 

where Vt is the (covariant) directional derivative, VtA^ = (p^VjA^ = A^ + r^jj^A'^ip^ [TflTi]. 
The gauge-invariant Mukhanov-Sasaki variables for the linearized perturbations, Q^ , obey 
an equation of motion with a mass-squared matrix given by [TJ [T3] 

M[j ^ g'' (VjVkV) - 7^^^M>^0^, (11) 

where Tl^ lmj is the Riemann tensor for the field-space manifold. 

We define adiabatic and isocurvature directions in the curved field space via the unit 
vectors a^ = 0^/o" and s^ = w^/w, where the turn-rate vector is given by uj^ = Vta^ , and 
oj = \uj^\. We also define slow-roll parameters [Tf [T^: 



H^ 



, f]^, = Ml ''''' ^^' \ r/,, = MpV ^ 



(12) 



Using Eq. (10), we have the exact relation, e = 3(t^/(o"^ + 2V). The adiabatic and isocur- 
vature perturbations may be parameterized as Tie = {.H/a)aiQ' and S = {H/&)sjQ^ , 
where TZc is the gauge-invariant curvature perturbation. Perturbations of pivot-scale /c^, = 
0.002 Mpc~^ first crossed outside the Hubble radius during infiation at time t*. In the long- 
wavelength limit, the evolution of TZc and S for t > t* is given by the transfer functions 

mm 



T-jis {t* , t) 



dt' 2uj{t')Tss{U,t'), 



tt 



Tss{t*,t) = exp 



dt' /3{t')H{t') 
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with 
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Given the form of T-jis, the perturbations TZc and S decouple if u^ = 0. 

The dimensionless power spectrum for the adiabatic perturbations is defined as V-jiik) = 
(27r)^^A;'^|7^cP and the spectral index is defined as n^ — 1 = dlnV-jz/dlnk. Around t^, the 
spectral index is given by [21 El CH [15] 

n,{Q = l-6e{Q + 2r]^^{Q. (15) 
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At late times and in the long- wavelength limit, the power spectrum becomes 

Vn = Vn{K)[l + Tls\ (16) 



and hence the spectral index may be affected by the transfer of power from isocurvature to 
iabatic modes: Usit) = rigit^) + H^^i^dT-jis/dt^) si 
The mass of the isocurvature perturbations is [7J 



adiabatic modes: n,(t) = n,(t,) + H^^dT-jis/dQ sin(2A), with cos A = T-jisil + T^^y^/^. 



,2 



^^l = 3H' (^r/,, + ^^ . (17) 

For fig < 3H/2 we have Vs{K) ~ V-jiiK) and hence Vs — 'Pii{t*)Tss at late times. 
In the Einstein frame the anisotropic pressure 11* • oc T'.- for i ^ j vanishes to linear 
order, so the tensor perturbations hij evolve just as in single-field models with Vt — 
128 [/J(t)/Mpi]2(A;/A;,)-2% and therefore [21 [H [15] 



^T 16e 

We demonstrated in [8J that for Higgs inflation, the system quickly settles into an ef- 
fectively single-field evolution independent of the fields' initial conditions, following a brief 
transient period. Here we find similar behavior even for potentials with nontrivial features 
like ridges and valleys: the system rapidly falls onto an effectively single-field, slow-roll 
attractor solution independent of the couplings A/ and ^i. 

We first consider the case in which the system inflates in a valley along the x = direction. 



perhaps after first rolling off a ridge. In the slow-roll limit and with X ~ X ~ 0, Eq. (10) 
reduces to [S] 



3v^5: 



Using Eq. ([8j) we may integrate Eq. (19), 

K2 4 
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A^*, (20) 



where N^, is the number of efolds before the end of inflation, and we have used 0(t*) 3> 



(pitcnd) — 0. Note that Eq. (20) is independent of A/. We arrive at comparable expressions 
if the system falls into a valley along some angle in field space, 9 = arctan {(f) /x)- 

1^=0 ^ QM^i^ 



Inflating in a valley along the x = direction, Eq. (9) becomes ct'^\x=o — 6Mpi0^/0^ upon 



using ^^ > 1. Using Eqs. (0, (|l9|), and (Q in Eq. (g we find 
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(=1) 



To estimate ?7o-o- we use the single-field approximation, rj„fj = 2e — e/{2He) = e — a /{Ha) + 



O(e^) [2J. If we first substitute Eq. (19) into the expression cr 1^=0 — 'oM^^cjr /(jr and then 



differentiate, we arrive at 

1 



All dependence on A/ and ^/ has dropped out of these expressions for e and ?7o-o- in Eqs. 



(21) and (22). For a broad range of initial field values and velocities — and independent of 
the couplings — this entire class of models should quickly relax into an attractor solution 
in which the fields evolve along an effectively single-field trajectory with vanishing turn- 
rate, u^ ~ 0. Within this attractor solution we find analytically e* = 2.08 x 10~^ and 
7]^^^ = -0.0165 for N^ = 60; and e, = 3.00 x 10"^ and r]^^^ = -0.0197 for iV, = 50. To test 
this attractor behavior, we performed numerical simulations with a sampling of couplings and 
initial conditions. We fixed A^ = 0.01 and ^(^ = 10^ and looped over A^ = {0.5,0.75, 1} A^, 
g = {0.5, 0.75, 1} X^, and ^^ = {0.8, 1, 1.2} ^0. These parameters gave a variety of potentials 
with combinations of ridges and valleys along different directions in field space. We set the 
initial amplitude of the fields to be v0o + Xo = 10 x max [^^ ,^x ] (i^ units of Mpi), 
which generically produced 70 or more efolds of inflation. We varied the initial angle in 
field space, Oq = arctan (0o/Xo); among the values 6*0 = {0,7r/6,7r/3, 7r/2}, and allowed for 
a relatively wide range of initial fields velocities: 4>o,Xo = { — 10 |0sr|)O, +10 |0sr|}) where 



()SR is given by Eq. (19). 



Typical trajectories are shown in Fig. [2^. In each case, the fields quickly rolled into a 
valley and, after a brief, transient period of oscillation, evolved along a straight trajectory in 
field-space for the remainder of inflation with u^ = 0. Across this entire range of couplings 



and initial conditions, the analytic expressions for e and r^o-o- in Eqs. (21)-(22) provide close 
agreement with the exact numerical simulations. See Fig. 2b. 

As we have confirmed numerically, trajectories in the single-field attractor solution gener- 
ically have cu^ ~ between t^, and tend (which we define as the time when e = 1, or d = 0); 
hence Ttis ~ for these trajectories. The spectral index Usit) therefore reduces to Usit^) of 
Eq. (15), and r reduces to r = 16e [1 + 0{T^g)] ~ 16e. Using Eqs. (21) and (22), we then 



find 
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and hence n, = 0.966 and r = 0.0033 for A'^* = 60; and n. = 0.959 and r = 0.0048 for 
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FIG. 2: Left: Field trajectories for different couplings and initial conditions (here for (j)o,Xo = 0). 

Open circles indicate fields' initial values. The parameters {X-^,g,^xi^o} are given by: 

{0.75 A^,A0, 1.2 e<^,7r/4} (red), {A^, A^, 0.8 ^0, vr/4} (blue), {A^,0.75 A<^, 0.8^^, 7r/6} (green), 

{A<^, 0.75 A(^,0.8^(^,7r/3} (black). Right: Numerical vs. analytic evaluation of the slow-roll 

parameters, e (numerical = blue, analytic = red) and ?]„„ (numerical = black, analytic = green), 

for A(^ = 0.01, A;^ = 0.75 A,^, g = \^, .^^ = 10^, and ^^ = 1.2 ^(^, with ^o = 7r/4 and 

4>o = Xi) = +1O|0sr|- 



A^* = 50. We also calculated n^ and r numerically for each of the trajectories described 
above, and found rig = 0.967 and r = 0.0031 for A^^, = 60, and rig = 0.960 and r = 0.0044 
for A^* = 50. These values sit right in the "sweet spot" of the latest observations. (See Fig. 
1 in jlj.) Even for a low reheat temperature, we find n^ within 2a of the Planck value for 
A^* > 38. The predicted value r ~ 10~'^ could be tested by upcoming CMB polarization 
experiments. 



For the running of the spectral index, a = dng/dlnk, we use Eq. (23), the general 



relationship (dx/dlnk) 1=,, ~ (i/H) \^, [2], and A^* = A^tot ~ /^ * Hdt to find 



to 
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which yields a = -5.83 x 10"^ for A^, = 60 and a = -8.48 x 10"^ for A^", = 50, fully 
consistent with the result from Planck, a = —0.0134 ± 0.0090, indicating no observable 
running of the spectral index P] . 

Meanwhile, for every trajectory in our large sample of numerical simulations we find 
|/nl| < 0.1, consistent with the latest observations [5]. In these models /nl is exponentially 
sensitive to the fields' initial conditions, requiring a fine-tuning of O(10~^) to produce |/nl| > 
1 [7]. In the absence of such fine-tuning these models generically predict |/nl| ^ C^(l)- 

Models in this class should produce entropy efficiently at the end of inflation, when 



Cii4>^Y < MpY and the couplings of the Jordan-frame potential dominate. For models with 
polynomial potentials of the form 0" near the minimum one expects Wcs = (w) = {n — 
2)/(n + 2) [1]. Eq. ([6| then suggests Wefi = 1/3. Tree-level masses may be added to V . For 
realistic values, mi <^ Mpi, the mass-terms do not alter the slow-roll attractor behavior [7], 
though they would dominate dynamics near the end of inflation and produce w^s = 0. More 
generally, inflation in these models ends with one or both fields oscillating quasi-periodically 
around the minimum of the potential. Thus preheating could be efficient |2], |16] . 

The models in this class predict three basic possibilities for isocurvature perturbations, 
depending on whether inflation occurs while the fields are in a valley, on top of a ridge, or 
in the symmetric potential of Higgs infiation. In the case of infiating in a valley, rjss > 1 
so nil H"^ > 9/4 and the (heavy) isocurvature modes are suppressed, T^s — > 0. In the case 
of infiating on top of a ridge, rjgg < so fJ^^/H"^ < and the isocurvature modes grow via 
tachyonic instability, Tss ^ 1- And in the case of Higgs infiation, fii/H^ — [8J, giving 
Tss ~ C(10~^). In each of those situations, meanwhile, o;^ ~ and T-jis ~ 0. Hence the 
primordial isocurvature fraction |3], 



B- (k)- ^'^^^ - ^^'^ f25l 



may easily distinguish between the various situations. In particular, if the system evolves 
primarily in a valley, then /3iso ~ for scales k corresponding to A^=k = 60 — 50. If the system 
evolves along the top of a ridge, /3iso ~ 1 across those same scales. Scenarios in which the 
fields begin on top of a ridge and roll off at intermediate times can give any value < /3iso < 1 
depending sensitively upon initial conditions. For Higgs inflation, Pi^o = 2.23 x 10~^ for 
N^ = 60 and Aso = 3.20 x 10-^ for N^ = 50. 

Multifleld models of inflation with nonminimal couplings possess a strong single-fleld at- 
tractor solution, such that they share common predictions for Ug, r, a, /nl, and for efficient 
entropy production across a broad range of couplings and initial conditions. The predicted 
spectral observables provide excellent agreement with the latest observations. These mod- 
els differ, however, in their predicted isocurvature perturbation spectra. We are presently 
conducting an expanded numerical study to understand the sensitivity of measures like /3iso 
on couplings and initial conditions. 
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